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We analyze a toy model which has a structure similar to that of the recently found QCD evolution 
equations, but without transverse dimensions. We develop two different but equivalent methods in 
order to compute the leading-order and next-to-leading order Pomeron loira diagrams. In addition to 
the leading-order result which has been derived within other toy models [3 j| 3|> we can also calculate 
the next-to-leading order contribution which provides the (agay) correction. We interpret this 
result and discuss its possible implications for the four-dimensional QCD evolution. 
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I. INTRODUCTION 

There have been major breakthroughs in understand¬ 
ing high-density QCD evolution recently. Several impor¬ 
tant observations have been made: (i) gluon number fluc¬ 
tuations have a big effect on the evolution towards gluon 
saturation Q, (ii) a connection between high-energy 
QCD evolution and statistical physics models of reaction- 
diffusion type was found 0, which has clarified the in¬ 
terpretation of fluctuations in an event-by-event picture, 
and (iii) it was realized that the relevant fluctuations 
are missed @ by the existing Balitsky-JIMWLK equa¬ 
tions 0,0- These observations have led to the creation of 
new QCD evolution equations which include gluon num¬ 
ber fluctuations or Pomeron loops 0, [ifl [HI O ■ 

The new equations describe the usual BFKL evolution, 
Pomeron splittings and Pomeron mergings and, thus, 
through iterations, Pomeron loops. A lot is known about 
the structure of the new equations mniiiiiiiiii, 
however, because of their complexity, little about their 
solutions. So far, analytical results for the energy depen¬ 
dence of the saturation momentum and for the scaling 
behavior of the T-matrix at asymptotic energies im and 
some numerical simulations of approximations to the ex¬ 
act equations [Hill are available. Solutions to the new 
equations in the range of collider energies are desired. 

In this paper, we consider a toy model which has a 
structure similar to that of the new QCD evolution equa¬ 
tions, but which has no transverse dimensions. This 
simplification allows an exact solution of the evolution 
equations. The particle-particle scattering amplitude is 
calculated. The effects of saturation and unitarity, or 
“Pomeron” loops which are the characteristics of the new 
evolution equations, are investigated. 

A decade ago, Mueller 0 studied a similar toy model 
without transverse dimensions which was motivated by 
the QCD dipole model ^3- He has been able to solve 
the toy model exactly and calculate the ^'-matrix. Ob- 
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servations made by studying this toy model (e.g. the 
one that fluctuations in particle numbers lead to a di¬ 
vergent multiple scattering series) have been shown nu¬ 
merically [^ 1^ to be valid also in the four-dimensional 
QCD. The hope that some of the results obtained in the 
toy model discussed in this paper may also apply in four¬ 
dimensional QCD is one of the motivations for this work. 

In the toy model proposed by Mueller 0, particle- 
particle scattering in the center-of-mass frame at total 
rapidity Y was considered. In this model, each parti¬ 
cle evolves to a dense system through particle branching 
(splitting). Unitarity effects (at leading order accuracy) 
are correctly described due to the multiple scatterings 
between the dense systems. There are no particle merg¬ 
ings in the wavefunctions of the particles in this model, 
therefore, the model fails to describe saturation effects 
in the wavefunctions of the particles which become im¬ 
portant at very high rapidities Y > Yg = -^ In In 
the language of loops, only the loops stretching over 
the rapidity interval greater than Y/2 (“large” loops) 
shown e.g. in Fig.QJI are included (at leading-order ac¬ 
curacy and for y < Yg). “Large” loops are created by 
matching the evolution (“Pomeron” splittings) of both 
particle’s wavefunctions at Y/2. However, there are no 
loops inside the particle’s wavefunctions (“small” loops) 
where the “Pomeron” s plit s and merges, as shown e.g. 
in Fig. [IP( see, e.g. ref.[^), within the rapidity inter¬ 
val Y/2. The above limitations of the toy model are, of 
course, properties of the QCD dipole model ^ as well 
and are well-known. 

Mueller’s toy model breaks down at U — Yg. In a later 
work Mueller and Salam 0 did include saturation effects 
by imposing boost invariance on the small-a; evolution in 
a toy model, getting an expression for the iF-matrix which 
is valid also for Y > Yg. The effect of saturation was an 
extension of the validity region of the result obtained for 
y < Ts in 0 also to the region where Y > Yg. 

Our toy model which is inspired by the new QCD evo¬ 
lution equations 0, [HIH [12 allows us to study unitar¬ 
ity effects as well as saturation effects in a more refined 
and systematic way. We calculate the leading order (LO) 
and the next-to-leading order (NLO) loop contributions 
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to the particle-particle scattering amplitude. The LO re¬ 
sult, which is the LO contribution of the “large” loops, 
agrees with the results from dll Our NLO result con¬ 
sisting of the NLO contribution of the “large” loops and 
the LO contribution of the “small” loops is new. The 
NLO correction with respect to the LO result is of order 
ttgOF. This is the main result of the paper. It tells us 
that the NLO contribution is negligible as compared with 
the LO contribution up to rapidities Y < l/aag. 

Our LO result agrees with the result derived within 
Mueller’s toy model up to T - Vs- Beyond Ys Mueller’s 
toy model fails to include saturation effects. The result 
in [2 which involves saturation effects is the same as our 
LO result until T ~ 1/ aa^ which is parametrically much 
larger than Yg. However, beyond Y ~ XjoLO^ the NLO 
corrections become important which are not taken into 
account in Mueller and Salam’s toy model 0. 

Presumably our toy model reveals the properties of the 
four-dimensional QCD evolution to some extent. In fact, 
its behavior is a heuristic sign that the NLO contribution 
can be neglected in the QCD evolution as well as long as 
Y <C llaal- Such a NLO correction seems to appear 
also in real QCD. Considering the one-loop diagram for 
example, the LO contribution in QCD is (a^) exp {2aY) 
which comes from varying the size of the loop from 0 to T, 
where a = ap — 1 = In 2. The NLO contribution 

(a^) aY exp (aY) comes from changing the location of 
a fixed zero-size loop from 0 to Y. The suppression of the 
NLO contribution by the factor aYje^"^ directly leads to 
the general correction a^aY (the same happens in our 
toy model). 

The paper is organized as follows: In Sec. QH we show 
the toy model. Two methods are developed in order 
to solve the toy model and some discussion is provided 
on the consequences of the NLO corrections. They are 
explained in Sec. inii In the Appendices we show the 
Lorentz-invariance of the toy model and give the solu¬ 
tion to another commonly used zero-dimensional model. 
Finally, we discuss the results and give the conclusions. 

II. A TOY MODEL WITH “POMERON” LOOPS 

In this section, we consider a toy model without trans¬ 
verse dimensions which has a structure similar to that 
of the recently found four-dimensional QCD evolution 
equation mill. The toy model describes the 

rapidity evolution of the particle number n of a system. 
The dynamics in zero-transverse dimensions is given by 
the following Langevin equation^ 

— = an — jJn^ + V2an v{y) (I) 


^ We provide the solution to the zero transverse dimensional 
sFKPP equation with a different noise term -y/2(Q!n — jSffl) i'{y) 
in Appendix. m 


where n(y) is a Gaussian white noise: {v{y)) = 0 and 
{v{y)v(y')) = S{y - y'). Eq. (HJ is known as the zero 
transverse dimensional “Reggeon field theory” equation 
and it contains the projectile-target duality (see Ap¬ 
pendix 0. The structure of Eq. o is obvious: the first 
term on the r.h.s. represents the growing of the particles, 
the second term describes particle recombinations which 
limit the growth at a maximum occupancy, and the third 
term describes the particle number fluctuations. These 
are essentially the main ingredients of the real QCD equa¬ 
tions. 

With Ito’s calculus, one can write Eq. o into an infi¬ 
nite hierarchy of coupled evolution equations, 

-= + k{k — , (2) 

dy 

where : ) is the expectation value of k- 

particles during the evolution and it should be under¬ 
stood as a normal ordered number operator (i.e. the 
factorial moment {n{n — 1) ■ ■ ■ {n — k + 1))) according to 
its physical interpretation. Here, kand^'> is the growth 
term, k{k—\)an'd^~^^ the fluctuation term, and 
the recombination term. Eqs. are the zero- 

transverse-dimensional analog of the real QCD equations 
(see A simple derivation and a discussion in the 

context of QCD evolution of eqs. d d can be found 
in . The solution to another widely discussed Langevin 
equation is shown in Appendix 151 


III. SOLUTION TO THE TOY MODEL 

Eq. d is quite complicated since is coupled with 
and in the evolution, i.e., Eq. d is just a 

particular equation within an infinity hierarchy of cou¬ 
pled equations. However, the following observation helps 
to solve Eq. @: one can neglect the recombination term 
if one starts with a dilute object at the beginning of rapid¬ 
ity evolution, on the other hand, one can see that = 
is an approximate fixed point of Eq. d in ex- 
tremely high rapidity and ^ << I limit. This will 

be more clear if one thinks about this in terms of the evo¬ 
lution equation of (one can define .): 

= a [A:T('=) - -F k{k - , in 

which it is easy to find limy^oo = I. Therefore, if 
one starts with a dilute object at Y = 0 (n <C = I/a^), 
then one can consider the evolution equation without the 
recombination term at the early stages of the evolution, 
and turn on the recombination term as a perturbation 
at the end of the evolution which forms the “Pomeron” 
loops. Following this philosophy, we have found two dif¬ 
ferent methods to solve Eq. @ which we will present in 
the following subsections. Before entering the calcula¬ 
tions, we should state that we always work in a paramet¬ 
rical limit in which —> 0 and aY oo, but 0 ^ 6 “^ is 
fixed and finite. 
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can be solved exactly if given some physically chosen ini¬ 
tial conditions. In this paper, we consider the evolution 
of a single particle for which the initial conditions are 


.(1) 


y=o 


= N and n. 


(k) 


k>l 


y=o 


= 0. The solution to Eq. @ 


with the above initial conditions reads 


nf \y) = IVfc!e'=“«(l - 


( 4 ) 


and corresponds to the A:-particle density obtained via 
splitting from a single particle after the evolution over 
the rapidity y. To turn on the perturbation, let us start 
with the first equation in the infinite hierarchy 


FIG. 1: Three kinds of reggeon graphs, in which the curly 
lines represent BFKL ladders. Diagram (A) is a simple BFKL 
ladder exchanged between the two particles (horizontal lines); 
Diagram (B) is the LO Pomeron loop diagram which has been 
calculated in ref. QilH; Diagram (C) is the NLO pomeron 
loop digram which was not included in the dipolar approach 
of ref. U12, 1^ while it can be computed in this toy model. 


A. Method I: Integral Representation 


Eq. @ without the recombination term, 


dn\ 


( 1 ) 


dy 


= 


( 5 ) 


which has the solution 


iW(y) = [ dye°‘^^ 

Jo 




aY 


1 


)• ( 6 ) 


dn. 


(fc) 


dy 


— = kan'^'^ + k{k — l)an^ , 


( 3 ) 


It is straight forward to recognize that the first term in 
Eq. ® is the usual BEKL term, and the second term cor¬ 
responds to the one-loop diagram with the right intrinsic 
minus sign. One can also see that is the NLO con¬ 
tribution in the parametrical limit we mentioned above. 


To get the LO two-loop diagram shown in Fig. ^B, we start with the second equation in Eq. 


dn\ 


( 2 ) 


dy 


= 2ari'^^ 2an''Q^ — 2Pn^\ 


( 7 ) 


whose solution is 


pY 

= n^^\Y)-2p dye 
Jo 




2Ne'^°‘^{l - e-“^) - 2 • 3!-Afe3“’^(l - ^ -b ...) , 


( 8 ) 


then turn on the perturbation, — Pri'i \ in order to take into account the recombination from to 

rpdd (in addition to the recombination from to given by Eq. O), which leads to 


.(1) 


pY 

(F) = n«(y)-/3 / dye 
Jo 


a{Y-y) „( 2 ) 


(y) 




3! ill 


fVe^“'^(I - 


4aY 

^aV 


+ ‘^(Pv))- 


( 9 ) 


The third term in Eq. represents the LO two-loop diagram shown in Fig. ^B. 

In order to calculate the NLO two-loop diagram shown in Fig. ^C, one has to enforce a recombination to happen 
after a splitting. Thus, starting with = iVe“^, one splitting given by ri' 2 '^ {Y) = 2a dye‘^°‘^^' followed 

by one merging = —/3 dye°‘^^~y^ n'^\y), yields the result for the one-loop diagram. A further splitting and 

a subsequent merging gives the result for the NLO two-loop diagram, 

= + , ( 10 ) 
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In fact, in order to calculate any loop diagram, one just needs to know the merging and splitting vertices which 
according to Eq. @ read, respectively, 

fY 

AnW(y) = -kf3 / , 

pY 

AnW(y) = k{k-l)a / 

Jo 

For particle-particle scattering, the (fc -I- l)-th order diagram contains fc-splittings and fc-mergings: The leading- 
order diagram corresponds to the diagram in which all fc-splittings occur before all the fc-mergings, the next-to-leading 
order diagram is always the diagram in which exactly one merging takes place before one splitting. The sum over all 
LO and NLO loop diagrams gives the LO and NLO contributions to the T-matrix, respectively. 

For particle-particle scattering, the sum over the LO loop diagrams (leading and next-to-leading contribution) reads 


( 11 ) 

( 12 ) 


'-LO 




l-{k-lf^+0 


while the sum over the NLO loop diagrams (leading contribution) is 


'-NLO 


(T) = ^(-l)''+^iV-fc!e'=“^ ^ 


(fc-l)(fc-2)2 av 


oaY 



(13) 


(14) 


The resulting ^-matrix, S = 1 — T (with T = which includes LO and NLO loop diagrams, can be 

written in the form 


OO 

S{Y) = E (-1)" • {l + alaY ■ {3k^ - fc) + h{k) + {aloYf ■ Mk) + ...} , (15) 

A;=0 


where we have put ^ in close analogy with the four-dimensional QCD and TV = 1. fi{k) and / 2 (fc) are functions 
of k, which are too complicated to be calculated in this method. The series in Eq. is a divergent series, however, 
it is Borel-summable. After the Borel-summation, the expression for the S'-matrix becomes rather complicated. 


SiY) = 


1 


■ exp 


-I- (alaY) 


■ exp 


r( 0, 
1 


1 


r 0 


4 + 


10 


+ {/i(ase“’^)} + (oi'^aYp |/2(a^e“^)| -I- ... , 


(a2e“i")2 


( a 2 e “^)2 


(16) 


however, it gets considerably simplified when a2e°‘^ becomes large 


S{Y) = 


1 

2 Q,y exp 


1 


r 0 


1 


l + 4,alaY + 0 


(17) 


where r(0, a;) is the incomplete gamma function. The NLO correction term 4a^aY inside the brackets with respect 
to the leading contribution is new. 


B. Method II: w-representation 

In this section, we develop another more powerful 
method, the ^-representation, to solve the toy model. 
This method exhibits the pole-structure of the toy model 
which bears some resemblance to the pole-structure 
of the BFKL equation in QCD. The ^-representation 
is equivalent to the integral representation method in 


Sec. 111! Al however, it enables us to go one step further by 
resumming all the a^aY terms in the S'-matrix. There¬ 
fore, this method allows us to get an expression for the 
S-matrix which is assumed to be valid up to the rapid¬ 
ity Yc 2 = where the next-to-next-to-leading order 

{p (a^) oT) becomes of order one. 

Let us show how to solve Eq. with the tu- 
representation method. One starts with the Laplace 
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transform of Eq. © 


a;n('=) (w) - (t) = (w) + fc(fc - (w) - fca^n('=+i) (w) , 




(18) 


where we have used the definitions t = ay and rS^\u}) = dtrS^\t)e With the initial condition n, 


( 1 ) 


= N 


and n, 


(fe) 


k>l 


y=o 


y=o 


= 0, Eq. 11811 can be written into a matrix forr 


(io-I ai 0 0 0 

-2 w - 2 2al 0 0 

0 ... ... ... 0 

0 0 —k(k — 1 ) (jj — k kal 

\ 0 0 0 . j 


V 


(i)(a;) \ 


(N\ 

(2) (a;) 


0 


= 

0 

W(u;) 


0 

... y 

V 0 / 


(19) 


It is fairly complicated to solve this equation directly because of the infinite dimension it has. Nevertheless, there are 
two ways to attack this problem. First, we can follow the idea that the is small, which reduces the importance of 
recombination at the beginning of the evolution. Then we can solve the above linear algebra equation perturbatively 
which will prove that this approach is completely equivalent to the integral representation method we developed in 
the previous section. Second, noticing that Eq. itTil has an exact solution at finite dimensions, we can truncate the 
Eq. 11911 to a k X k linear algebra equation and then take the k oo limit to get the solution of our infinite dimension 
matrix. 

First, let us turn off the perturbation and only solve the unperturbed matrix to get the corresponding nQ^^(w). 
One can easily find the solution to be y_ k) ~ ^ After the inverse Laplace 

transformation, 


r(<^) 


(fe)/ ^ 1 


^S+icxD 




( 20 ) 

( 21 ) 


the result is exactly the same as what we have obtained in the previous section. 

Turning on the perturbation, we find that the first-order correction (w) to (w) satisfies the equation: 


/w-1 0 0 0 

-2 w - 2 0 0 

0 ... ... 0 0 

0 0 —k{k — 1) u — k 0 

V 0 0 0 . ) 


0 \ ( AnJ^^(w) \ 

0 


AnJ^^(a;) 
An^^^(a;) 


/ _„,2„(2) 


/ 


^s"’0 


(w) \ 




-ka^,n‘^o^^\uj) 


( 22 ) 


The general solution of this equation is too complicated 
to be written down. However, the strategy of how to get 
it is straightforward. For example, we can easily obtain 

AnJ^^(a;) = which corresponds to the one- 

loop diagram. It yields the same result as given in Eq. 
after the inverse Laplace transformation. This example 
shows that the result for each diagram calculated in the 
previous section becomes pretty much simplified in the 
^-representation. To see this more clearly, one can do the 


Laplace transform of Eq. CD and Eq. 03 which yields 
— ka'^ 

An«(a;) = -f An('=+i)(a;) , (23) 

ijj — k 

AnW(w) = . (24) 

According to these equations, one can directly write down 
the contribution for any diagram in this w-representation 
since the merging of A: -|- 1-particles to fc-particles is noth¬ 
ing but adding a pole with a prefactor —kal in 
ta-complex plane and the splitting of fc — 1-particles to 
A:-particles is nothing but adding a pole with a pref- 
























actor k {k — 1). Although the problem has been simpli- Eq. 11911 is nevertheless unavailable this way. 
fled a lot in the ^-representation, the exact solution to 


6 


We get the solution to Eq. ca as follows: We first truncate the infinite dimension matrix m into a finite k x k 
matrix, then exactly solve the amputated equation. Let us show the solutions for the first four fc-values: 


k = I case: One obtains which transforms to nl^\aY) = . 

k = 2 case: One obtains — (cj- i){ij- 2 )+ 2 a‘^ which transforms to 

n\^\aY) = . 

^ ^ ^ 1 - 4a2 1 - 4a? 


(25) 


7 0 /-\ 1 • (1) / \ A^[(a; —2)(cc;—3) + 12a^l . . . 

• k = 6 case: One obtains = ?—ttt—— Q^, o 2 /? —?vr- inen one can set the determinant 

oX3V / (cj —l)(u; —2)(u;—3)+2 q;|(7u;—9) 

(w — 1) (w — 2) (w — 3) -I- 2a^ {7lu — 9) = 0 which yields three solutions. Each of those three solutions stands for 
three poles uji = l-|-2a^, uj 2 = 2-|- lOa^ and W 3 = 3 — 12ag, respectively. The inverse Laplace transform changes 
the integer exponents to non-integer renormalized exponents. 




1 + 3Q:^ (l+2a?lay _ 2 a^ -f lOOal f2+10a?)ay , 72af (3-12a?lay 
1 -I- a? 1 — 14a? 1 — 31a? 


(26) 


• fc = 4 case: One obtains The denominator has 4 poles 

which are wi = 1 -I- 2a? -I- o (a?), uj 2 = 2 + 10a? -I- o (a?), W 3 = 3 -I- 24a? -I- 1164a? and W 4 = 4 — 36a? — 1080a?. 
Then, the inverse Laplace transform gives 


nW(ay) = iV(l + o(a?))e^+2“')“'^-2!iVa?(l + o(a?))e(™'K 

YSlNai (1 -h o (a?)) - 7344 (a?) (27) 


From the above calculation for n(^)(ay), one can already 
see that for an arbitrary k, the first k — \ poles are stabi¬ 
lized to the corresponding solutions ojj = jY (3j^ — j)a?-I- 

O (a?) (1 < j < k) with a fixed coefficient j\N {oilY ^ 
and the last pole is always LOk = k—k{k — 1 )^ a!?-|-C> (a?) 
with a wrong coefficient since it misses the recombination 
from [k + 1 ) particle mode as a result of truncation at 
the k particle mode. In addition, it is easy to see that 
the (3j^ — j)a? term is corresponding to the sum of the 

I 


NLO contributions (j-I-1) j^a?ay of the leading loop 
diagrams (e.g. Fig. ^ B) and the leading contributions 
(—j (j — 1)^ a?ay) of the NLO loop diagrams (Fig.^C). 
On the other hand, the —k {k — 1)^ a? term in the last 
pole only corresponds to the leading order contributions 
of the NLO loop diagrams (e.g. Fig. Q] C) as a conse¬ 
quence of the truncation at dimension k which rules out 
the leading loop diagram but allows the NLO loop dia¬ 
grams at the k particle mode. 


The pattern of the exact solution is transparent from the above discussion and the generalization to infinite di¬ 
mensions follows by simply taking the limit k 00 . The exact solution for the S'-matrix in infinite dimensions 
reads 


Se.act{aY) = (-1)'^ («")" (l + O («?)) . (28) 

This result agrees with what we have found in the previous section, see Eq. (CHI, however, it goes a step further since 
the k (3fc — 1) a?ay term is now in the exponent which can be understood as a resummation of all a^aY terms via 
exponentiation in Eq. 03- _ 


We notice that the S'-matrix is not Borel-summable since is beyond the Borel-summability accord- 
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ing to Nevanlinna’s theorem . However, one can gen¬ 
eralize the Borel-summation technique and show that the 
above series still has a finite sum by using analytic con¬ 
tinuation. Generally speaking, one can see that the S- 
matrix is the analytic function of two independent vari¬ 
ables, aj. and aY, when is negative in the complex 


plane. Since = 0 is not a singularity, one can use 
analytic continuation to define the other half plane when 
aj. becomes positive. In the following, we provide a gen¬ 
eralized Borel-summation method which yields a finite 
result. 


We use a technique together with a subsequent Borel-summation in order to maneuver the divergent series Eq. 
into a definite and well-defined sum. Let us define u = and 7 = Sa^aY, then we obtain 


S'(aF) = 


00 


= ^(—1)^ J dbb^e (y\f^ J dxexp + 2jxk)^ , 

r-\-oc ^ noc / h\ 

= \ — / dx exp (—70;^) — / ddexp (-) 

V J-00 uJo \ uj 


1 -h be+'^'y^ ’ 


^2 d, exp (- 7 x>) ^ exp r (^0. 


1 


yg-e 27 x 


(29) 

(30) 

(31) 

(32) 


To check the convergence of this integral, one can expand T (0, z) at z ^ 00 which gives T (0, z) ~ 
showing that the above integral is definitely finite and well-defined. The x-integration can not be performed analyti¬ 
cally. However, one can use the saddle point approximation to evaluate the above x-integral, in which one finds the 
saddle point x = — 1 , and 

g,2g(l-7a2)ay) ^ (O’ g^2g(l-7al)ay 


v2p(l-4a2)ay 


ue 


+2-yx 


(33) 

(34) 


A good agreement of this result with the numerical eval¬ 
uation of Eq. m is shown in Eig. ifM . 

Eq. 121 is a fairly exact result in the high-energy limit 
since it presumably resums all the sub-leading diagrams 
up to o?^olY level, only neglects O (a^) (in the prefactor) 
and O (a^aY) ( in the exponent) corrections. If we com¬ 
pare this with Eq. d of the previous section, we can see 
that Eq. d is an approximate solution which results 
from Eq. m when a^^oY <C 1. The w-representation 
is equivalent to the integral representation method in 
Sec. IIII Al however, it enables us to go one step further 
by resumming all the a^^oY terms in the S'-matrix. Thus, 
the S'-matrix in Eq. m is assumed to be valid up to the 
rapidity Yc 2 = where the next-to-next-to-leading or¬ 
der (O (agOy)) becomes of order one. 


C. The Consequences of the NLO Correction 

In Mueller’s toy model Q, which was inspired by the 
QCD dipole model (DM) S'-matrix for particle- 

particle scattering in the center-of-mass frame is given 


by 

gDM^Y) = f2 (I) (I) > (35) 

where Pn(Y) is the probability density for having n par¬ 
ticles in the wavefunction of an initial single particle af¬ 
ter the evolution up to Y. Eor large rapidities, Pn{Y) 
reads Q 

J^(n,E) = ^ (36) 

n{Y) 

with n{Y) = exp(asF). Using u = nln{Yj2) and v = 
m/n{Yf2) and changing the summations in eq. II35II into 
integrals (this change is rigorous since du=dv= l/n(r/ 2 ) 
is infinitesimal) for large rapidities, one obtains 

/ OO nOO 

du / dv exp [—(E) uv — u — u](37) 

V(V/ 2 ) -(y/ 2 ) 

which after the integration over u and v leads to 






■In 




(1 -I- 


(38) 























FIG. 2: The S'-matrix as a function of aY for = 0.02 
in logarithmic scale: The solid line represents the solntion 
of the zero-transverse dimensional Kovchegov equation. The 
dashed line corresponds to the leading order result, the dotted 
line stands for the saddle point approximation and the 

dot-dashed line is extracted from the numerical evaluation of 
Eq. Idlll . 



FIG. 3: The T-matrix as a function of aY for = 0.02: The 
solid line represents the solution of the zero-transverse dimen¬ 
sional Kovchegov equation. The dashed line corresponds to 
the leading order result, the dotted line stands for the saddle 
point approximation lEU and the dot-dashed line is extracted 
from the numerical evaluation of Eq. lISTl . 


This result reduces for F <C Ks = — In ( ^ ) to 

<^3 \ / 

( 39 ) 

while for F ^ Fg it gets 


SDM^Y) 


1 


In 



(40) 


It is easy to see that the result in Eq. agrees with 
our LO contribution given in Eq. m only when F < Fg. 
For Y Yg the result from Mueller’s toy model given 
in Eq. inni) can not be trusted since it does not include 
saturation effects within the interacting particle wave 
functions which start becoming important at F ~ Fg. 
Mueller’s toy model fails to reproduce the LO contribu¬ 
tion of the LO diagrams (e.g. Fig^B) at F Fg, it 
does not include the NLO contribution of the LO dia¬ 
grams (terms including (fc— l)^aF/e“^ in Eq. 1131) 1 and 
it also misses the NLO loop diagrams (e.g. FigHJC). 

Mueller and Salam have been able to calculate the 
probability distribution P„(F) which includes saturation 
effects by imposing boost invariance on the small-x evolu¬ 
tion. They calculated the boost invariant S matrix in the 
rest frame of the target, and obtained the same LO result 
in Eq. nrii . We also hnd that it is interesting to calculate 
the S matrix in the center of mass frame. Inserting their 
result, 


-P*(«: y) = exp 


a^n — ay 


at V 


- 1 


^41) 


in Eq. (1^ . converting the summations into integrals, 

changing the integration variables from m — 1 to m and 

2 

from n — \ to n, setting e"® = 1, and using the integration 
variables v = njn{Yl2') and u = mln{Y 12), one ends up 
with the boost-invariant result for the ^'-matrix in the 
center of mass frame. 


^^(F) 


poo poo 

/ du dv exp [—Ogn (F) MW 
Jo Jo 

11 1 

_ - _no_ - _1 


.ln[a^e“n. 


u — ^2) 

(43) 

(44) 


This result agrees with our LO contribution in Eq. CZJ. 
Note, however, that the Mueller-Salam toy model does 
not incorporate the NLO correction given in Eq. HU. 
Our toy model, Eq. 03, tells us that the NLO correction 
is negligible for Y < Y^ = This means that the 

result in Eq. m is valid over a much larger rapidity 
range, up to F ~ F^, = as compared to the result in 
JSHl which is inspired from the dipole model and which 
breaks down at Fg = ^ In Therefore, the effect 

of saturation is to replace the result in Eq. (HHll by the 
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result in Eq. even when Y > Yg. It is interesting 
to note that the saturation effects are mathematically 
included by setting the lower limits in Eq. 123 to 0 as 
can be seen by comparing Eq. m with Eq. without 
modifying the probability distribution. (In addition, one 
can easily prove that the S matrix of dipole approach is 
boost invariant if and only if the lower limits in Eq. m 
are set to 0.) 


IV. CONCLUSION 

We have investigated a stochastic toy model without 
transverse dimensions (equivalently, a infinite dimension 
hierarchy of evolution equation) which naturally gener¬ 
ates Pomeron loops. We have computed the Pomeron 
loop diagrams to the NLO using two different methods. 
We find that the LO calculation agrees with the exisiting 
results in the literature P, 01 ■ The study of NLO graphs 
in this toy model generates the [a^aY) corrections. 

Observations made by investigating Mueller’s toy 
model 0 have been shown numerically by Salam 
to be valid also in the four-dimensional QCD. Pre¬ 
sumably our toy model reveals properties of the four¬ 
dimensional QCD evolution to some extent. We be¬ 
lieve that the NLO correction in the four-dimensional 
QCD would be of order o^aY as well (here we define 
a = ap — 1 = In 2 to be the BFKL pomeron inter¬ 
cept). The behavior of our toy model is heuristically in¬ 
dicating that NLO contributions can be neglected in the 
four-dimensional QCD evolution as long as V ^ W as 
well. This believe is strengthened even more by noticing 
that the NLO correction a^aY seem to naturally appear 
in QCD: There are two types of NLO contributions. 

• The first type of NLO contribution comes from 
NLO contribution of LO graphs(see example 
Fig.njB). For the one-loop diagram, for instance, 

the LO contribution in QCD is (a^)^ exp (20:1^) 
which comes from varying the size of the loop from 
0 to Y. The NLO contribution (cis)^ aY exp {aY) 
comes from changing the location of a fixed zero- 
size loop from 0 to Y. The NLO contribution is sup¬ 
pressed by the factor aYje°‘"^ with respect to the 
LO contribution, as in the toy model (see Eq. m), 
leading to the general form a^aY. 

• The second type of NLO contribution comes 
from LO contribution of NLO graphs(see example 
Fig.dC). For the two-loop diagram Fig. ^C, for 
instance, the LO contribution (ag)^aV exp (20:1^) 
comes from shifting the location of the fixed zero- 
size connecting Pomeron between two loops from 0 
to Y while the sizes of these two loops are changed 
in the same time. Since the total amplitude of this 
graph does not depend on the location of the con¬ 
necting Pomeron, one can obtain an extra factor 




FIG. 4: Graphs for one object scattering with three objects. 
Guriy lines represent BFKL ladders. In graph (A) a sin¬ 
gle particle evolves through Pomeron splittings and scatters 
off three particles, while in graph (B) three particles evolve 
through Pomeron mergings and scatter off a single particle. 


of aY by integrating from 0 to F. The NLO con¬ 
tribution (see Eq. CD) is suppressed by the fac¬ 
tor aYle°‘'^ with respect to the LO contribution of 
Fig.^B, as in the toy model (see Eq. C21)), leading 
to the general form o^aY. 

Therefore, we can see that the NLO correction gener¬ 
ally takes the form of oi\aY. A QCD model in which all 
Pomeron splittings occur before all Pomeron mergings in 
a boost-invariant manner (as in this toy model) would 
lead to the right LO result which would be valid up to 
the limit 
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APPENDIX A: THE EQUIVALENCE OF 
SPLITTING AND RECOMBINATION 

In this part, we show that Eq. 13 gives Lorentz in¬ 
variant results. To show this, let us start with a sin¬ 
gle particle, then split it to three particles which yields 
n^\y) = jV3!e^“^(l -- e““^)^ according to Eq. 0). Fol¬ 
lowing Kovchegov |^, the scattering amplitude is 

k^l 

in which (0) = ( and the ^ can be understood 
as the fact that those fc-particles in the target are the 
same. Therefore, it is very easy to see that As{aY) = 
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N (1 —which corresponds to the diagram 

shown in Fig.^jA. 

On the other hand, we can put the evolution into 
the target, let them merge into one object before scat¬ 
tering with the projectile (Fig. 2]B). This means that 


one starts with an initial condition n. 


(3) 


TO 


V=0 


y=o 


= Ns and 


= 0 in solving the Eq. After merging 

twice from n^\y) = one obtains = 


Np 




— e Finally, it is straight forward 


to see that Ajn{aY) = Ns (1 ~ e after set¬ 

ting ^ = Og. Presumably, Ag{aY) should be same as 
Ajn{aY) since they actually describe the same process, 
and they are indeed the same once we set N = Ns- 


Moreover, it is tempting to replace the target by a 
large nucleus with A nucleons inside it. Then, we need to 

define tW(0) = (a , and indeed one can see that 
As{aY) = Ajn{c(Y) as long as we set Ns = N (A^'j . 


I 

APPENDIX B: SOLUTION TO A SLIGHTLY DIFFERENT STOCHASTIC DIFFERENTIAL EQUATION 


We notice that there has been a lot of theoretical discussion an d numerical simulation on a slightly different 
stochastic differential equation which has a noise term y^2 (an — namely the stochastic Fisher-Kolmogorov- 

Petrovsky-Piscounov(sFKPP) equation. For completeness, we provide the solution to this equation by using the 
w—representation we developed above. The dynamics of this equation is given by two random processes: by an 
increase dY in rapidity, a particle can split into two particles with some rate a (A AA A -|- A) or two particles can 
merge into one with a rate 2/3 (A -|- A A) . From the equation. 


— = an — + \J2 (an — j3n^)v{y), (Bl) 


one can get 


—;— = (ka — k{k — 1) (3) n^^^ + k{k — (B2) 

dy 

Similarly, one can write Eq. (El into a matrix in the w-representation, then get the similar S'-matrix of the scattering 
between two dipoles. 


S{aY) 


^ (-1)'= fc!e(fc+2fc^aJ)aY+o(aj)aF ^ ^ ^^2^^ 


A;=0 


^2g(l-2a2)„y exp (^2g(l-4a?)aY) ^ ' 


r( 0 , 

1 


,2p(l-4a2)ay / ’ 


(B3) 

(B4) 

(B5) 


where S = 2alaY and v = The difference of between this result and the one in Eq. (IMIl is at the NLO 

correction, however, one can see that Eq. o has a more transparent physical meaning in QCD since it contains the 
projectile-target duality and one can easily relate it to the z ero transverse dimension reggeon theory. _ 
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